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Abstract 

New data on the J ftp suppression both in proton-nucleus and in lead- 
lead interactions have been presented recently by the NA50 collaboration. 
We show that these data, together with the final ones on sulfur-uranium 
interactions, can be described in the framework of the comovers interaction 
model with a unique set of three parameters : the nuclear absorption cross- 
section, the comovers interaction cross-section and a single (rescaled) absolute 
normalization. Expectations for J/ip suppression at RHIC are also discussed. 
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1 Introduction 



Before the Quark Matter conference of 2002, the NA50 interpretation of the 
data on J /if} suppression was as follows [1-3]. The pA, SU and peripheral Pb Pb 
data (up to Et ~ 354-40 GeV) can be described with nuclear absorption alone, with 
an absorptive cross-section a a bs = 6.4 ± 0.8 mb. At Et ~ 40 GeV there is a sudden 
onset of anomalous suppression, followed by a steady fall off at larger E T . However, 
at variance with this view, the most peripheral points in Pb Pb collisions lied above 
the NA50 nuclear absorption curve - which extrapolates pA and SU data. 

Two important sets of new data have been presented recently j3] The new 
NA50 data on pA reactions at 450 GeV/ c indicate a smaller value of a a b s than the one 
given above. However, within errors, pA and SU data can still be described with 
a single value of the absorptive cross-section a a b s = 4.4 ± 0.5 mb - substantially 
lower than the previous one 0. The new, preliminary, Pb Pb data [3], taken in 
2000 with a target under vacuum, are consistent with previous ones except for 
the most peripheral ones - which are now lower and consistent with the nuclear 
absorption curve In this way, the NA50 interpretation remains valid. However, 
the new data lend support to the interpretation based on comovers interaction - 
according to which some anomalous suppression is already present in SU collisions. 
Indeed, in the comovers approach 2 the sudden onset of anomalous suppression due 
to deconfmement is replaced by a smooth anomalous suppression due to comovers 
interaction. The effect of the comovers turns out to be negligibly small in pA but it 
is sizable in SU interactions. With the smaller value of cr a 5 S from the new pA data, 
there is more room for comovers in SU. 

The purpose of this work is to study the consistency of the new data on pA and 
Pb Pb interactions, together with the final SU ones, with the comovers interaction 
model jH] We proceed as follows. Since the effect of the comovers suppression is 
sizable in SU, but negligibly small in pA, a a bs has to be determined from the pA data 
alone. In previous works [8-9] we have used a value a a b s = 4.5 mb, as a compromise 
between NA38/NA51 [TO] and E537 [TT] pA data. Actually, it has been shown in [T2| 
that the old pA data are also consistent with <j a b s = 4.5 mb. As mentioned above, 
2 For reviews on deconfining and comover interaction models see UJ. Alternative models have 
also been proposed 7 . 
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this value has now been confirmed by the recent NA50 data 0, and will be used 
throughout this paper. The second parameter in the model, the comovers interaction 
cross-section a co , can then be determined from the centrality dependence of the J /if) 
suppression in Pb Pb collisions. Obviously its value is correlated with that of the 
third parameter of the model, the absolute normalization. This normalization, in 
turn, is strictly related to the one in SU. The ratio of the Pb Pb to the SU 
normalizations is equal to 1.051 ± 0.026 pj- This is a rescaling factor which takes 
into account both isospin and energy corrections. In this way, with a a b s fixed, the 
model is strongly constrained. 

The main drawback of the comovers model [8-9] was precisely a mismatch of 
about 30 % between the absolute normalizations in SU and Pb Pb [T2|. The origin 
of this mismatch is the following. The high values of the most peripheral Pb Pb 
data in the former NA50 analysis required a value <j co = 1 mb. As stated above, in 
the new data the most peripheral Pb Pb points are substantially lower and require a 
lower value of a co = 0.65 mb. Indeed, a smaller value of a co (with <r abs fixed) leads to 
a flatter centrality dependence of the J/if> suppression. This change in a co induces 
a change in the absolute normalization - which is now in good agreement with the 
(rescaled) one obtained in SU. 

The plan of this paper is as follows. In Section 2 we present a short summary of 
the comovers interaction model [8-9]. In Section 3 we apply it to Pb Pb collisions, 
where the data allow an accurate determination of both the comovers cross-section 
and the absolute normalization. We also compute the correlation between and 
Ezdc ~~ the energy of the zero degree calorimeter and discuss the J /if) suppression 
in the Ezdc analysis. In Section 4 we show that the pp, pA and SU data can be 
described using the same values of a a b s and a co as in Pb Pb and a single (rescaled) 
normalization - obtained from either SU or Pb Pb data. Section 5 contains our 
conclusions and expectations for J /if) suppression at RHIC. 

2 Comovers interaction in the dual parton model 

The cross-section of minimum bias {MB), lepton pair (DY) and J /if) event 
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samples are given by 

l£? B (b) oc o AB (b) (1) 

oc Jd 2 sa AB (b)n(b,s) (2) 

1 ab ( b ) / cPs <tab(&) n(6,s) S ahs (b,s)S co (b,s) . (3) 



Here ctab(&) = {1 - exp[-cr pp AB T AB (b)]} where T AB (b) = J d 2 sT A (s) T B (b - s), 
and T A (b) are profile functions obtained from the Woods-Saxon nuclear densities 
[T3j . Upon integration over b we obtain the AB total cross-section, cab- n(b,s) is 
given by 

n(&, s) = AB a pp T A (s) T B (b - s)/a AB (b) . (4) 

Upon integration over s we obtain the average number of binary collisions n(b) = 
AB a pp T AB (b)/a AB (b). 

The factors S a b s and S co in (j^J are the survival probabilities of the J/ip due to 
nuclear absorption and comovers interaction, respectively. They are given by [S] [H] 

Q abs (h x _ [1 - exp(-AT A (s) aabsW ~ exp{-B T B (b - s) a abs )} 

{ , ) ~ <s AB T A (s) T B (b - s) (5) 



S 00 ^, s) = exp 



o , co , , i ■-> - , 

Vdt 



CO 



-^NZ T (b,s)£n^^ 



(6) 



In (jBJ), Ny° DT (b, s) is the density of charged comovers (positives and negatives) in the 
rapidity region of the dimuon trigger and Nf = (3/irR 2 > )(dN/dy) y *=o = 1.15 fm~ 2 
[HJiniE] is the corresponding density in pp. The factor 3/2 in © takes care of the 
neutrals. In the numerical calculations we use a abs = 4.5 mb. The value of a co and 
the absolute normalization will be determined from the data. 

In order to compute the density of comovers we use the DPM formalism de- 
scribed in [03 ■ It turns out that the density of charged particles is given by a linear 
superposition of the density of participants and the density of binary collisions with 
coefficients calculable in DPM. All details can be found in jH] and [To] . 

Eqs. (0) to © allow to compute the impact parameter distributions of the 
MB, DY and J/ip event samples. Experimental results are plotted as a function of 
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observable quantities such as Ex - the energy of neutrals deposited in the calorimer. 
Using the proportionality between Et and multiplicity, we have 

Mb) = I Q N-Jb) . (7) 

Here the multiplicity of comovers is determined in the rapidity region of the Et 
calorimeter. The factor 1/2 is introduced because N co is the charged multiplicity 
whereas Et refers to neutrals. In this way q is close to the average transverse 
energy per particle, but it also depends on the calibration of the calorimeter. The 
correlation Et — b is parametrized in the form [3J E] 

P(E T , b) = . 1 exp {-{Et - E T (b)) 2 /2qaE T (b)} . (8) 

^2irqaE T (b) 1 J 

The Et distributions of MB, DY and J /if) are then obtained by folding Eqs. ©-© 
with P(E T ,b), i.e. 

I^ DY ^{E T ) = j d% I^ DY ' J/ \b) P(E T , b) . (9) 

The parameters q and a are obtained from a fit of the Et distribution of the MB 
event sample 3 . Note that since N™ (b) is nearly proportional to the number of 
participants (see Fig. 1 of jH]), our fit is practically identical to the one obtained 
[Hj using the wounded nucleon model. Actually, we obtain identical curves to the 
ones in Fig. 1 of ref. jUj - where the Et distributions of MB events of 1996 and 
1998 are compared with each other. The values of the parameters for the 1996 data 
are q = 0.62 GeV and a = 0.825. For the 1998 data, the tail of the Et distribution 
is steeper, and we get q = 0.62 GeV and a = 0.60 4 . In the following we shall use 
3 Note that the same value of the parameter a is used in the MB, DY and J/ip event sample. 
A priori there could be some differences in the fluctuations for hard and soft processes. Actually, 
it has been claimed in Ref. [H] that there is a small shift in Et between minimum bias, on one 
hand, and J/ip or Drell-Yan pair production on the other hand - induced by the dimuon trigger. 
However, this is of no relevance for the present work, since, so far, the only 2000 data available 
have been obtained in the so-called standard analysis - in which the genuine ratio of J/i/j and DY 

cross-sections is measured. 

4 At first sight these sets of values look very different from the ones used by the NA50 collab- 
oration. Nevertheless, they reproduce the same Et distribution. This is due to the fact that the 
product qa, which according to Eq. (JHJ, determines the width of the distribution, is very similar 
in the two cases. As for the difference in the values of q it is just due to its definition, which is 
different in the two approaches (Eq. Q, in our case). 
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the latter values. Indeed, according to the NA50 collaboration [2], the 1996 data 
(thick target) at large Et are contaminated by rescattering effects - and only the 
1998 data should be used beyond the knee. 

The model described above allows to compute the Et distribution of MB, DY 
and J /if} event samples between peripheral AB collisions and the knee of the Et 
distribution. Beyond it, most models, based on either deconfinement or comovers 
interaction, give a ratio of J/if> to DY cross-sections which is practically constant 
- in disagreement with NA50 data. A possible way out was suggested in 0. The 
idea is that, since Et increases beyond the knee due to fluctuations, one can expect 
that this is also the case for the density of comovers. Since does not contain 

this fluctuation, it has been proposed in to introduce the following replacement 
in Eq. © : 

N ZA b > s ) - 0> s ) = N ZA b > s ) F ( b ) ( 10 ) 

where F(b) = Et j 'Ej<(b) . Here Et is the measured value of the transverse energy and 
Et(P) is its average value given by Eq. (JJJ) - which does not contain the fluctuations. 

3 J/i/j suppression in Pb Pb 
a) Et analysis 

The new data jl] for the ratio of J/if> over DY cross-sections versus the energy 
of the Et calorimeter are shown in Fig. 1. They are compared with the results of 
the comovers interaction model described in Section 2. As explained there, there 
are two free parameters in the model (<r a fc s = 4.5 mb has been fixed) : the comovers 
interaction cross-section a co (which controls the centrality dependence of the ratio) 
and the absolute normalization. A good description of the data is obtained using 
<j co = 0.65 mb and an absolute normalization of 47. 

The only difference between our result and the one in || resides in the value of 
a co . Since the effect of the comovers increases with centrality, a larger (smaller) value 
of a co leads to a larger (smaller) variation of the ratio of J / if) over DY cross-sections 
between peripheral and central collisions. As mentioned in the Introduction, in the 
new NA50 analysis the values of this ratio for peripheral collisions are smaller. In 
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order to describe the new data, the value of a co has to be reduced. The curve in 
Fig. 1 corresponds to a reduction of a co from 1 mb (used in 8J) to 0.65 mb. 

Since the values of a co and of the absolute normalization are correlated, the 
decrease of a co induces, in turn, a decrease of the absolute normalization. While 
in jH] the value of the absolute normalization was about 30 % higher [12] than in 
SU, the one in Fig. 1 is in good agreement with the SU one. This will be shown in 
the next section. It is interesting that almost the same value of a co [a co = 0.7 mb) 
was obtained in from an analysis of SU data and former Pb Pb data pQ which 
covered a much smaller centrality range. In jTH] the absolute normalizations in SU 
and Pb Pb were in good agreement with each other. 

The DY cross-section in Fig. 1 has been integrated in the dimuon mass range 
2.9 to 4.5 GeV. Since the J/ip peak is inside this range, a model is needed in 
order to determine the DY cross-section. In the SU analysis, the GRV parton 
distribution functions at leading order (LO) have been used. Therefore, in order to 
use the (rescaled) absolute normalization of the SU data in Pb Pb (or vice versa), 
the same GRVLO distributions have to be used in the latter. This is the case in 
Fig. 1. In Pb Pb collisions, NA50 has also analyzed the data using, instead, MRS 
43 distributions. They have found 4J that in this case the absolute normalization 
is lower by about 10 %. The comparison of the comovers model with the data is 
presented in Fig. 2. The absolute normalization is 43. The values of <r a bs and a co 
are, of course, unchanged. 

b) E ZDC analysis 

The NA50 collaboration has also measured the J/ip suppression in Pb Pb as a 
function of the energy of the zero degree calorimeter (Ezdc)- In Fig. 1, the results of 
this analysis have been plotted as a function of Et, using the measured correlation 
between average values of E T and E Z dc- We see from Fig. 1 that the data obtained 
in the two analysis are consistent with each other, even for very central events, 
beyond the knee of the Et and Ezdc distributions. This important result has been 
predicted in [T2j. Its physical origin is the following. 
The energy of the zero degree calorimeter is given by 

E Z Dc{b) = [A — n A {b)]E in + a n A {b) E in . (11) 
7 



Here n^(6) is the average number of participants at fixed b : 




(12) 



A — n,A(b) is the number of spectator nucleons of A and Ei n = 158 GeV is the beam 
energy. While the first term in the r.h.s. of Eq. (fTTj) gives the bulk of Ezdc, the 
latter corresponds to the contamination by secondaries emitted very forward 
- assumed to be proportional to the number of participants, riA{b). Here also the 
value of a can be precisely determined from the position of the "knee" of the Ezdc 
distribution of the MB event sample measured by NA50 We obtain a = 0.076 



Eqs. ((Zj) and (fTTj) give the relation between b and Et and b and Ezdc, respec- 
tively. These relations refer to average values and do not contain any information 
about the tails of the E T or E Z dc distributions. Eqs. (J7|) and (fTTj) also lead to a cor- 
relation between (average values of) Et and Ezdc- This correlation (12] gives a good 
description of the experimental one [3] . It is practically a straight line 5 and therefore 
can be accurately extrapolated beyond the knee of the E T and E Z dc distributions. 
It turns out that this extrapolation describes the measured E T — E Z dc correlation 
quite well 6 - even for values of Et and Ezdc in the tails of the distributions. This 
result suggests a correlation between b and Ezdc of the form 



Folding (0) and (fTBf) we obtain the E Z dc distribution of MB events. It describes 
fT2"j the one measured by NA50, not only up to the knee, but also in the tail of the 
distribution. This result shows that the J /ip supression versus Ezdc is just obtained 
from the corresponding one versus Et by applying the Ezdc — Et correlation, 
even for very central events beyond the knee of the distributions. In the 1996 and 
1998 NA50 data, the J/ip suppression versus Ezdc indicated some features ("snake 
shape") not present in the ones versus Et- Such differences are no longer present 
in the new data. 

5 This is due to the fact that Ny° a (b) in Eq. J7J) is practically proportional to 71,4(6) (see Fig. 1 
of®). 

6 One can understand the physical origin of this extrapolation if one assumes that a fluctuation 
in Et is essentially due to a fluctuation in ua - which, in turn, produces a corresponding fluctuation 
in Ezdc, via Eq. (fTTfl . 



P(Ezdc, b) = P(E T , b)5 (E T - E t (E Z dc) - 



(13) 
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4 J/ip suppression in pA and SU 

Let us compute next the ratio R of J /if} over DY cross-sections in SU at 
200 GeV/c per nucleon in our model. We use, of course, the same values of the 
parameters as in Pb Pb : a a bs = 4.5 mb and a co = 0.65 mb. To get this ratio versus 
b, the only new ingredient is the multiplicity of comovers - which is again computed 
in DPM, in the way described in Section 2. To compute R(Et), we also need the 
Et — b correlation in SU, which is parametrized as in Eq. (jSJ). The parameters q 
and a have been obtained from a fit of the /^-distribution of DY given in ^H] • We 
obtain q = 0.69 GeV and a = 1.6. (R(Ezdc) has not been measured in SU). In 
577, data do not extend beyond the knee of the i^-distribution. Therefore, effects 
such as E T fluctuations, Eq. (fTUj). are not relevant here. 

Our results are shown in Fig. 3. We see that the Ex dependence of the sup- 
pression is reproduced. This indicates that there is, indeed, room in SU for the 
(comparatively small) suppression by comovers. As discussed above this could also 
be inferred from the different central values of a a b. s obtained in pA and SU. The 
absolute normalization of the curve in Fig. 3 is 45. Thus the normalizations in Pb 
Pb and SU are consistent with each other. This normalization is 4 % smaller than 
the one obtained from the Pb Pb data - in perfect agreement with the rescaling 
factor discussed in the Introduction. 

In pA collisions, the effect of the comovers is negligible and, therefore, the de- 
scription of the pA data is the same as in the NA50 analysis [5], since, as discussed 
above, the value of a a b s they obtain is practically identical to ours. The correspond- 
ing normalization is 20 % higher than the one we have obtained in SU. This is also 
consistent with the expected rescaling factor between the two systems - which takes 
into account the difference in energy as well as in the rapidity regions covered by 
the dimuon trigger. 

5 Conclusions and outlook 

The NA50 deconfining scenario has been described in the Introduction. In this 
work we have presented a different scenario in which the sudden onset of anoma- 
lous suppression due to deconfinement is replaced by a smooth one resulting from 
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comovers interaction. This anomalous suppression is already present in SU and 
peripheral Pb Pb collisions. 

We have presented a description of the NA38/NA50 data on the J/ip suppression 
in pp, pA, SU and Pb Pb interactions in a comovers model. The model is strongly 
constrained by the existing data. This can be seen in the following way. In SU, the 
effect of the comovers is rather small and, thus, once the value of a a bs = 4.5 mb is 
fixed, the absolute normalization depends little on the exact value of a co . Since the 
normalizations in SU and Pb Pb are strictly related, we are left with a single free 
parameter, a co , to determine the J/ijj suppression in Pb Pb (both in absolute value 
and centrality dependence). The model is, thus, strongly constrained and provides 
a unified description of the data in the various systems. On the other hand, it is 
not possible to describe the former NA50 data on Pb Pb collisions for peripheral 
events in a consistent way. Indeed, these data require a value a co = 1 mb. As 
shown in ^2] , this, in turn, leads to a mismatch of about 30 % between the absolute 
normalizations in SU and Pb Pb systems. Furthermore, with equal normalizations 
in SU and in Pb Pb, the J / ip suppression is always larger in the latter than in the 
former, even for very peripheral events (see Fig. 7 of [12J). 

Let us discuss briefly the expectations for J/ip suppression at RHIC in the co- 
movers interaction model. The calculation of the survival probability S co is quite 
safe. Indeed, since a co is a cross-section near threshold, the same value obtained at 
SPS should be used at RHIC. The situation is quite different for S a b s - Many au- 
thors assume that a a b s is the same at RHIC and at SPS. It has also been suggested 
that it can be significantly larger at RHIC. However, it seems plausible that at mid- 
rapidities, nuclear absorption at RHIC is small due to the fact that, contrary to SPS, 
the cc pair is produced outside the colliding nuclei. It is therefore crucial to have 
data on J/ ip production in pA interactions at RHIC. If S a b s ~ 1, the J / ip suppression 
at RHIC and SPS will be comparable, since the smallness of the nuclear absorption 
will be approximately compensated by the increase of the comovers suppression - 
due to a larger comovers density at RHIC. Very preliminary data tend to indicate 
that this is indeed the case. Detailed calculations will be presented elsewhere. 

A quantitative analysis of the new NA50 data in the deconfining scenario is still 
missing. On the other hand, the centrality dependence of the average pt of J/ip is 
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better described in the comovers approach than in a deconfining scenario j^. At 
RHIC energies, a small nuclear absorption mpA collisions (i.e. S a b s ~ 1), would be a 
very interesting situation in order to discriminate between the comovers interaction 
model and a deconfining scenario. Indeed, in the latter, the shape of the centrality 
dependence would be almost flat for peripheral collisions (below the deconfining 
threshold) and would decrease above the threshold. Such a behavior would be a 
clear signal of deconfinement. On the contrary, in the comovers scenario, the fall-off 
would be continuous, from peripheral to central collisions, and determined by the 
same value of a co obtained from CERN SPS data. 
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Figure Captions 



Fig. 1 : Ratio of J/ip to DY cross-sections versus E? in Pb Pb collisions at 158 
GeV/c per nucleon (solid line). The preliminary data are from jlj. GRVLO parton 
distribution functions have been used in order to calculate the DY cross-section in 
the mass range 2.9 to 4.5 GeV. 

Fig. 2 : Same as in Fig. 1, using MRS 43 parton distribution functions. 

Fig. 3 : The ratio of J / ip to DY cross-sections as a function of E? in SU collisions 
at 200 GeV/c per nucleon. The data are from pCHj ■ 
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